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Abstract. We count ]B, C]-graiiied, fc-factor integers which are simuhaneously 
-B-rough and C-smooth and have a fixed number k of prime factors. Our aim 
is to exploit expHcit versions of the prime number theorem as much as possible 
to get good explicit bounds for the count of such integers. This analysis was 
inspired by certain inner procedures in the general number field sieve. The 
result should at least provide some insight in what happens there. 
We estimate the given count in terms of some recursively defined functions. 
Since they are still difficult to handle, only another approximation step reveals 
their orders. 

Finally, we use the obtained bounds to perform numerical experiments that 
show how good the desired count can be approximated for the parameters of 
the general number field sieve in the mentioned inspiring application. 

Keywords. Smooth numbers, rough numbers, counting, prime number theo- 
rem, general number field sieve, RSA. 
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Let us call an integer A-grained iff all its prime factors are in the set A. Then 
an integer is C-smooth iff it is ]0, C]-grained, and B-rough iff it is oo[- grained. 
You may want to call an integer family coarse-grained if it is ]B, C]-grained and the 
number of factors is bounded. We consider the number of integers up to a real positive 
bound X that are ]B, C]-grained and have a given number k of factors, in formulae: 

(0.1) A^ci^) ■■=#{n<x\3pi,...,pk £¥n]B,C]: n = pi---pkY 

We always assume that C > B, since otherwise the set under consideration is empty. 

Such numbers occur for example in an intermediate step in the number field sieve 
when trying to factor large numbers. During the sieving integers are constructed as 
random values of carefully chosen small-degree polynomials and made i?-rough by 
dividing out all smaller factors. The remaining number is fed into an elliptic curve 
factoring algorithm with a time limit that should allow to find factors up to C. To tune 
the overall algorithm it is vital to know the probability that the remaining number is 
C-smooth. Assuming that the polynomials output true random numbers the counting 
task we deal with is the missing link since estimates for counting i?-smooth numbers 
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are known, see for example the overview by iGranvilld (l2008l ) . It is difficult, however, 
to compare our results to existing ones, since we consider integers with a fixed number 
of factors. This restriction is completely sufficient for the application we have in mind. 
To our knowledge these kind of investigations cannot be found in the literature. Our 
main result implies the following 

Corollary. Fix k G N>2, a > and e > 0. Then for large B and C = B^~^^ we 
have uniformly for x G [B^{1 + e), C'^(l — e)] that 

Actually, we can describe a piecewise smooth function that approximates vr^ up to 

an additive error of or der O (^'^^ tlis hidden constant depends on e, k and a 
only, see [Theorem 9.21 [Theorem 9.31 

To avoid difficulties with non-squarefree numbers, instead of numbers we count 
lists of primes 

(0.2) K%^c (^) := # {(Pi, • • • ,Pfc) G (Pn ]5, C])*^ \pi---pk<x]. 

It turns out that there are anyways only a few non-squarefree numbers counted 
by vr^f;, namely kl ■ ir'^ ^ (x) ~ k^^{x). This would actually be an equality if 
there were no non-squarefree numbers in the count. We defer a precise treatment 
until [Section 8l 

We head for determining precise bounds ^ (x) or vr^ ^, (x) that can be used in 
practical situations. However, to understand these bounds we additionally consider 
the asymptotical behaviours. This is tricky since we have to deal with the three 
parameters B, C and x simultaneously. To guide us in considering different asymp- 
totics, we usually write B = x^ ^ C = x"' and 7 = /3(1 -|- q). In particular, C = 
So we replace {B,C,x) with (x,/3,7) or (x,a,/3), similar to the considerations when 
counting smooth or rough numbers in the literature. Alternatively, it seems also 
natural to fix x somehow in the interval ]B^, C'^] by introducing a parameter ^ by 

X = B'^^^C^ = 

Now the parameters are {B,C,£,) or {B,a,£,). 

As a first observation, note that ^ (x) is constantly for x < B^ and constantly 
(7r(C) — ir{B))^ for x > and grows monotonically when x goes through [B'',C^]. 
Here vr denotes the prime counting function. For 'middle' x-values the asymptotics 
can be derived from our main result: 

Corollary. Fix k G N>2, a > and e > 0. Then for large B and C = we 
have uniformly for x G [B''{1 + e), C'^(l — e)] that 



<c.W€e(j^)=e(^) 
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Later, we specify a piecewise smooth function ^ which approximates with 

an error of or der O (^'^^ w here the hidden constants depend on e, k and a only, see 
[Theorem 3.21 [Theorem 3.51 and subsequent considerations. 

The preceeding resuh is in contrast to the asymptotics at x = C^: 



This observation is explained as follows: Note that roughly half of the numbers up 
to X are in the interval [\x, x\ and similarly for primes. Thus the behaviour of those 

candidates largely rule '^_b,c X- For an x = B^-^C^ with a fixed ^ G ]0,A;[, it 
is mostly determined by the requirement that the counted numbers are i?-rough, 
and we thus observe a comparatively large fraction of C] -grained numbers. In 
the extreme = C^, most candidates are ruled out by the requirement to be 

C-smooth and thus we see a much smaller fraction of ]B, C]-grained numbers. 

The case that the intervals ]B^, C^] are disjoint for considered values k is especially 
nice, as then the number of prime factors of a S-smooth and C-rough number n can 
be derived from the number n. So we assume in the entire paper that C < B'^ for 
some fixed s > 1. (Clearly, we cannot have a fixed s that grants disjointness for all 
intervals. But for the first few we can.) In the inspiring number field sieve application 
we have C < iJi-232Tji] 'pj^jg ensures that th e intervals ] B^ , C * ^] for k < 5 are disjoint. 

A further application is related to RSA. Decker Sz Moree ( 20081 ) give estimates for 



the number of RSA-integers. They use one ad-hoc definition proposed by B. de Weger. 
However, it is not so clear which numbers we should call RSA-integers. A discussion 
and further calculations to adapt our results to the differe nt shape are needed. We 



have treated these issues in iLoebenberger &: Niiskenl (j201ll ). 



As our basic field of interest is cryptography and there the largest occurring 
numbers are actually small in the number theorist's view, we assume the Riemann 
hypothesis throughout the entire paper. We use the following version of the prime 
number theorem: 



Prime number theorem 0.3 (jVon Kochlll90ll . ISchoenfeldlll976l ). If (and only if) 
the Riemann hypothesis holds then for x > 2657 



|7r(x) — li(x)| < —^/xlnx, 
Svr 



where li(x) = dt. 



^Side remark: to indicate how a real number was rounded we append a special symbol. Examples: 
TT = 3.14J, = 3.142V = 3.1416T = 3.14159JI. The height of the platform shows the size of the left- 
out part and the direction of the antenna indicates whether actual value is larger or smaller than 
displayed. We write, say, e = 2.72T = 2.7in as if the shorthand were exact. 
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We ha ve numerically verified this inequality for x < 2^*^ ~ 1.1 • 10^^ based on lKulshal 's 
tables (IKuls ha 2008) and extensions built using the segmented siever implemented by 



Qliveira e Silvai (|2003i) . We are confident that we can extend this verification much 
further. In the inspiring application we have x < 2^'' and so for those x we can take 
this theorem for granted even if the Riemann hypothesis should not hold. 
We arrive at the following description of the desired count: 

Theorem. Let B < C = -8^+" with a > and fix k > 2. Then for any (small) 
e > and B tending to infinity we have for x G [B^{1 + e),C^{l — e)] a value 



a G 



k\{l+a)'' ' k\ 



^B,C 



with Cfc 
ix)-2 



mini 2-4^, 2 



-k 



(fc-1)! 



X 



InB 



< (2'' 



such that 

X 



B B 



□ 



Also without assuming the Riemann hypothesis we can achieve meaningful re- 
sults provided we use a good unconditional version of the prime number theorem 
wh ose error estimate is at least in O (j;^^ ■ The famous work bv Rosser &: Schoen- 
feld (|l962l . I1975I ) is not sufficient. Yet, iDusartI (|1998l ) provides an explicit error 
bound of order O fj^^^j a-nd Ford ( 2002a ) provides explicit error bounds of order 

though this only applies for x beyond 10^^^ or even much 



O (^x exp 



A(lnz)3/5 
' (lnlna;)i/5 



later depending on A and the O-constant, see lFact 6.11 

1. The recursion 

The essential basis for the analysis of the counting functions ^, is the following 
simple description. 



Lemma 1.1. For all k £ N>n we have the recursion 



E 



^B,C 



PfcePn]B,C] 



based on 



^B,c 



(x) 



Proof. In case k > we have 



(x) 




Jo if x£[0,l[, 
\l ifxe [l,oo[. 

n]B,c])''\pi---pk<x] 

Pk = p 



Pk& 
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The case /c = is immediate from the definition. 



□ 



From the definition (0.2) or from lLemma l.ll it is clear that 



ifxG[0,5[, 
7r{x)-TT{B) iixe[B,C[, 
[tt{C)-tt{B) ifxe[C,oo[. 



This reveals that the case distinction in (j leaves its traces on higher kJ^C- -^^'^ 
further calculations it is vital that we make this precise. This will enable us later to 
do our estimates. For /c S N we distinguish k + 2 cases: 



X is in case (/c, —1) :<;=^ x G 
x is in case {k,j) :<;=^ x S 
X is in case (/c, h) :<;=^ x € 



0,5' 



where j E N<fc. Note that most cases are characterized by the exponent of C at the 
left end of the interval. 



Lemma 1.2. For k,j£N,0<j<k and x in case {k,j), we have 



E 



4,c i^/Pk) + Yl 



"3,0 



where in the first sum x /pk is in case {k — l,j — l) and in the second in case {k — l,j). 
For j = —1, and j = k we do not split the sum, as then all x/pk are in one case 
anyways. For j = the left part is zero, so that the splitting there is less visible. 

Proof. We only have to verify that x/pk G [B'^^^^^C^ ,B''^^^'^C^^^[ for pk G 
Fn]B,x/Qk-i-j(jj] and x G [B'^'iC^ , B''-'^-^C^+'^[. Similarly, the statement for 
the second sum is established. □ 

For example, we obtain 

fO if X G [0,^2 [, 

Ep,em]B,^] Ep^gpn]B,-|] ^ ^ ^ [B^BC[, 
x) = { Ep2ePn]f ,c] I^p,6Pn]B,^] ^ 
+ Z]p2ePn]B,f ] X^piePn]B,C] ^ 

.J2p2&r\]B,c] J2pi&n]B,c] 1 if X G [C2,oo[. 



(1.3) 



'^B,C 



if X G [BC, [, 
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So for (J (x) we have four cases with four 2-fold sums. In general (j (x) has 
k + 2 cases with 2^ /t-fold sums. Well, we better stop unfolding here. 

Based on the intuition that a randomly selected integer n is prime with proba- 
bility jji^ we can replace X^pgpnj^ c] fip) ^ith dp. This directly leads to the 
approximation function. We prefer however to follow a better founded way to them 
which will also give information about the error term. 

2. Using estimates 

From the recursion for ^ (x) it is clear that we have to compute terms like 

pePn]B,C] pePn]B,C] 

To get good estimates for such a sum we follow the classic path, as Rosser &: Schoen- 
feld (1962): we rewrite the sum as a Lebesque-Stieltjes-integral over the prime count- 
ing function tt{p). Then we substitute 7r(t) = Li(i) -|- E{t), keeping in mind that we 
know good bounds on the error term E{t) by the lPrime number theorem 0.31 Finally, 
we integrate by parts, estimate and integrate by parts back: 

p£¥n]B,C] 

fit) dnit) 

B 

C fC 
B 

C fu\ rC 



f{t)dU{t)+[ fit)dE{t) 
JB 

Is ^ ^t+fiC)E{C) - f{B)E{B) - 1^ Eit) dfit) . 



The existence of all integrals follow from the existence of the first. If the sum kernel 
/ is differentiable with respect to t we can rewrite E{t) df{t) = f'{t)E{t) dt. 
Now we can use the estimate on the error term E{t): 



p&r\]B,c] 

< \f{C)\E{C) + \f{B)\E{B) + r \f'{t)\E{t) dt. 

JB 

What remains is, given the concrete /, to determine the occurring integrals. For 
the counting functions ^ — as one would guess — this task is more and more 
complicated the larger k is. Clearly, smoothness properties of / must be considered 
carefully. 

During all this we make sure that the involved functions stay sufficiently smooth: 
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Lemma 2.1 (Prime sum approximation). Let f, f, f be functions R>o 
that f and f are piecewise continuous, f + f is increasing, and 



s.>o 



such 



f{x)-f{x) <f{x) 



for X G M>o. Further, let E{p) be a positive valued, increasing, smooth function of p 
bounding \tt{p) — Li(p)| on [B,C]. (For example, under the Riemann hypothesis we 
can take E{p) = ^^Inp provided p > 2657.) Then for x G M>o 



peFn]B,C] 



X) 



where 



m = [ 

Je 

9{x) 



fjx/p) 

B 

^ f{x/p) 



dp, 



\np '^^^'^^^ + f)i^/B)E{B) + (/ + /) {x/p)E'{p) dp. 



c 



Moreover, g and g are piecewise continuous, and g + g is increasing. 



Proof. The assumption immediately implies 



pePn]s,C] pePn]s,C] 



< 



pem]B,c\ 



Using the techniques just sketched we obtain 'Yl,pevr\]B C\ f^''^ Ip) ~ Ib ^^\Jp^ dp + 
Jb fi^/p) d-E'(p) with E{p) = Tr{p) — Li(p). Shifting the second term to the corre- 
spondingly transformed error bound, we now have 



E f(^/p)-£ 



pem]B,c] 



f{x/p) 
Inp 



dp 



< 



MA Ap+ (/ + f){x/p) dE{p) 



This bound always holds, yet we still have to estimate E{p) in it. Abbreviating 
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h{p) = (/ + f){x/p) we estimate the last integral: 

h{p) AE{p) = h{C)E{C) - h{B)E{B) - [ E{p) dh{p) 

B J B 

< h{C)E{C) + h{B)E{B) - [ E{p) dh{p) 

Jb 

= + h{C)E{C) + h{B)E{B) 
- h{C)E{C) + h{B)E{B) 

rC 



+ [ h{p) dE{p) 
Jb 



For the inequality we use that h(p) is decreasing in p. Collecting gives the claim. 

Finally, g and g being obviously piecewise continuous it remains to show that 
^ + ^ is increasing. By the (defining) equation 

{9 + 9)ix) = 1^ (7+ f){x/p) + E'ip)^ dp+2(7+ f){x/B)E{B) 

this follows from f + f and E being increasing. □ 

3. Approximations 

Based on ILemma 2.11 we recursively define approximation functions ^ and error 
bounding functions k^,^- Recursive application of ILemma 2.11 will lead to a good 
estimate in [Theorem 3.21 For the understanding we further need to determine the 
asymptotic order of the functions defined here, which we start in the remainder of this 
section. [Theorem 3.51 relates the functions ^ and to some easier manageable 
functions. These in turn are computed or estimated, respectively, in [Section ^ and 
[Section 5l 

Definition 3.1. For x > we define 

^B,c (x) := Kb,c (x) , ^B,c (x) ■■= 
and recursively for k > 

+ 2(75^ + ^^) {x/B)E{B) 



+ («B,(7 + ^Bfi) {x/pk)E'{pk) dpk . 
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These functions now describe the behavior of ^, nicely: 
Theorem 3.2. Given x G M>o and k e'H. Then the inequality 

4,c {x)-liB,c \ <Kb,c (.x) 

holds. 

Proof. Using ILemma 2.11 the claim together with the fact that (j + q is 
increasing follow simultaneously by induction on k based on c ~ ^% C ^'^'^ '^B C ~ 
0. ' ' ' □ 

In order to give a first impression we calculate k]^ q and k}^ q . Analogous to 
ILemma 1.21 we split the integration at / B^~^~^C^ that the parts fall entirely 
into case (A; — 1, j — 1) or into case {k — 1, j): 



7 _B'=-i-JC-' 



Also for Q this can be done, simply split the occurring integrals at / B^^^^^ ■ 
Now unfolding the recursive definition of k;]^ ^ gives: 



ifxG[0,5[, 

(^) = <! /I isk d^^i ifxG[i?,C7[, 
Jb ink dpi -lixdiCcol 

{0 ifxG[0,S[, 
E{x) + E{B) ifxE[i?,C[, 
E{C) + E{B) ifxe[C,oo[, 

which corresponds exactly to the approximation of the prime counting function by 
the logarithmic integral Li. In case k = 2 we obtain 



(3.3) 







fB p 

JB JE 



1 



IB Inpi lnp2 

rC rP2 1 



dpi dp2 
dpi dp2 



JB In Pi In P2 

+ /J Ib 

Jb Ib lnpi\np2 dpi 



if X e [0,^2 [, 

if X E [B^,BC[, 

if X E [5C, C2 [, 

if X E [(7^,00 [. 



This is now exactly the transformed version of (1.3) as announced there. You may 
ask where you can find a display of the error term corresponding to (3.3) Well, we 
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have computed it. But the resulting terms are so complex that we didn't really learn 
much from it. Here is an expression for x G [B^,5C[: 

x_ 

Jb ^"^P 

Though this term is still handleable, it becomes apparent that things get more and 
more comphcated with increasing k. We escape from this issue by loosening the bonds 
and weakening our bounds slightly. The first aim will be to obtain easily computable 
terms while retaining the asymptotic orders, the second aim will be to still retain 
meanmg ful bounds for the fixed values B = 1100 • 10^, C = 2^7 - l, e {2, 3, 4} from 
our inspiring application. 

Our next task is to describe the orders of (j and (j. The main problem in an 
exact calculation is that most of the time we cannot elementary integrate a function 
with a logarithm occurring in the denominator. But using B < p < C we can obtain 
a suitably good approximation instead by replacing with in the integrals. At 
this point we start using C < B^ and rewrite C = where a is a new parameter 

(bounded by s — 1). For the time being you can consider a as a constant, but actually 
we make no assumption on it. This leads to the following 

Definition 3.4. For x > we let 

X''{x):=K%^c{^), A°(x):=0, 
and recursively for k > 

+ 2 (a^-^ + A'^-i) (x/B) E{B) 

+ r (a^-1 + A*^-^) {x/pk)E'{pk) dpk. 



We observe that InS < Inp^ < InC = (1 + a)lnS. Thus we obtain ^ d^?G 
Ib positive integrable function /. By induction on k we 



i+a' 
obtain 
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Theorem 3.5. Write C = and fix k ^ N>o. Then for x G M>o we have 

1 



4,C {x) G 



^B,c [x) e 



1 



1 



A'^ {x) , 

A^' (x) . □ 



In order to determine at least the asymptotic orders of ^ and /t^ ^, (along with 

precise estimates) it remains to solve the recursions for and A'^, or at least to 
estimate these functions. As a first step, we rewrite all integrals in IDefinition 3.41 in 
terms of g defined by pk = B^~^^°'. 

IDefinition 3.41 continued. Rewrite x = B^^^'^ with a new parameter ^ and let 
}k Jk (^fc+?a^ ^^^\k ._ \k (Actually, you may think off (0 := 

jk^j^k+^a-^ for any family of functions f^.) 

Lemma 3.6. For k > {) we now have the recursion 



(0 = a A'^-i (e - g) 5^+^" d^ , 

^0 



A' (0 = « 



Jo 

+ 2 (a'^-i + A'^-i) {OE{B) 

+ a\nB (\^-^ + A'^-i) - g) ^'(5i+f")5i+e« d^. . 



□ 



4. Solving the recursion for A 



k 



To construct a useful description of the function A'^ we make a small excursion and 
consider the following family of piecewise polynomial functions. 

Definition 4.1 (Polynomial hills). Initially, define the integral operator M by 

iMfm= C f{i-g)dg= fig)dg 
Jo J^-i 

for any integrable function f : M — )• M. ATow for k G N>i let the k-t\i polynomial hill 
be 

m'' := Mm''^^ 

based on the rectangular function fh} given by m}{£^) = 1 for ^ G [0, 1[ and m^(^) = 
otherwise. 
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Actually, = Mm? if we let mP = 5 he the 'left lopsided' Dirac delta distribution 
defined by its integral ^ 5{t) dt being for ^ < and 1 for ^ > 0. Contrastingly, 

the standard Dirac delta function is balanced and has 6{t) dt = ^. However, we 
will stick to the lopsided variant throughout the entire paper. By T)^ we denote the 
differential operator with respect to ^ . 

Lemma 4.2 (Polynomial hills). (i) m''{^) = for ^ < or ^ > k. 

(Hi) fn!' restricted to [j,j + 1[ is a polynomial function of degree k — 1 for j G 
{0, . . . ,k — 1}. In particular, is smooth for ^ S M \ {0, . . . ,k}. 

(iv) fh^ is {k — 2)-fold continuously differentiable. 



(v) Conversely, the conditions (i) through (iv) uniquely determine m . 



(vi) P^m'^ = AiV^^m^ ^ as long as i < k — 2 and even for i = k — 1 when read for 
distributions. 

(vii) The function m'^ is symmetric to |; rh^(^) = m^{k — ^). 
(via) For ^ G + 1[ (corresponding to case {k,j) ) we have 

(ix) The next derivate can only he correctly described as a linear combination of 
Dirac delta distributions: 



0<j<k 

(x) For any < i < k we obtain the following explicit description: 



(xi) Further, m^{i) = -^m^'^iC) + ^m^-^{i - 1) holds. 



Curiosity: The function ^/k-m^ can also be described as the volume of a slice through 
a (k — l)-dimensional unit hypercube of thickness orthogonal to a main diagonal. 

That's the same as saying it is times the (k — l)-volume of the cut between a 
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if 


1 
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-5) 


-1 


5 ^ 





1 2 3 



6g^-60g-^ + 210g^ -300€ + 155 




3(4-0 + 4 



4(5-5) -5 



1 2 3 4 ■ 
-1 




Figure 4.1: Graphs of the polynomial hills for k < 5 and their derivatives 



fc-dimensional hypercube and the hyperplane X]i<j<fc Qi = This last interpretation 
makes it obvious that 

Jo 

since this is the volume of the /c-hypercube. 



Proof. From the definition (i) through (vii) follow directly. Further, note that 
{MV^f) (0 = fiO - /(e - !)• This is obvious from (Mfm = /(q) dg. 

We prove (viii) by induction on k. For k = 1 this is true by definition. So consider 
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k > 1. The recursion for m differentiated k — 1 times yields for G + 1[ 



V 



k-2~k-l, 



V. 



k-2~k-li 



{-ly 



k-2 



i-1 



, YA: - 1 



(ix) follows similarly. 

To prove [(^consider h{0 = ^^^^.^^^j " 0'"'- Then Pj-^/i = 

V^f^mf' using [(vm)| And obviously P|/i(0) = = 2?|m'=(0) for < ^ < A; - 1, so 
that inductively (with falling I) we get = V^in^. 



As we do not need (v) and (xi) we leave these proofs to the interested reader. □ 
Most of the following is easier if we first renormalize A'^' . So we let 



(4.3) 



1 



The recursion for A now turns into 



'\k _ X/ix^-^ 

'^norm ^ ^norm ■ 



Theorem 4.4 (Approximation order). For any ^ G M we have 
Jo Jo 



1 



{-alnB) 



"norm 



o<j<UJ 



(-1) 



0<^<A;-1 

j cutexpj^. {—{^ — i)aln B) 



i-alnBY 



where cutexp^(C) = exp(^) — J2o<e<k-i W ~ Tl,i>k T\ - '^^^ ^^P^^ss cutexp^ 
using the incomplete Gamma function T{k,() = e"^ g'^^^ dg by cutexp^(C) = 

exp(-C)r{A;,0) • 
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Proof. The definition for A'^ turns into 

Jo 

Now, since Ai commutes with this integration this immediately imphes that 

^ /"OO 

Jo 

which is the first stated equahty noting that m*^ is zero outside [0,k]. By partial 
integration we obtain 

amB amB Jq 

^^^(alnS)' ^ {alnBfJo 



Using the description of P|m^ from Lemma 4.2(ix) the last integral turns into 

the claimed sum of the second stated equality. Expressing D^~^~^m^{^ — g) us- 
ing Lemma 4.2(x) and rearranging slightly yields the third equality. □ 

We are going to estimate the estimation of the error in the next section. To 
that aim we first need to estimate X^q^.^- If /c > then, based on [Theorem 4.41 and 
fh^iC) < 1, we obtain the upper bound 



~ f'°° 1 

ALm(0= / B-^^m''{^-g)dg< — 
Jo « In 



a 111 B 

For A; = we have A°o,ni (0 = for ^ > and so X'^^^^ < 1 will do for all 

We first describe the qualitative behaviour of X'^^j.^. Actually its graph looks like 
a slighty bias wise hill. 

Lemma 4.5. The function X^^^.^ (^) = B'^" B^°'rh''{g) dg is zero at ^ = 0, posi- 
tive at S, = k, more precisely 

~ /I — B^°^ ^ ^ 

Anorm (^) = (^^]^ 

and there is a position ^\ G ]0, A;] such that it is increasing on ]0, [ and decreasing 
on Further, Cl>\. 

2 2 
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Proof. First, inspecting 



Jo 



if C < 0, 



alnB 
1-B 




if^>l. 1 

shows that for k = 1 all claims hold with = 1. So in the remainder of this proof 

2 

we assume k > 1. 

Next, compute Aji|orm (^) inductively: 

Jo Jo 

= / 5^^*" da •-B-(^-i)" / B^'^m'^-^T) dr = 
Jo Jo 



1 - B-^ 

a In S 



Here we have substituted t = g — qj. and collapsed the new integration interval 
[—gk,k — gk\ for r to [0, k — \] since m^~^ vanishes on the difference Moreover, based 
on A^„rm (0 = B-^"^ BS'^m^ig) dg we obtain 



(4.6) P^Ai^^ (0 = -a\nB- X^,,^ (0 + m\i), 



and infer that T>^\^^^^ {k) = — alni? ^ alnB j negative. 
Finally, compute the derivate of A^orm differently 

T^Arm iO = ^5 (^^ B-^^m^C -g)dg^= B'^^ j^^ B^^V^m\g) dg . 

The integral kernel B^^T>^rh^{g) is positive on ]0, |[ and negative on Thus 
increases on ]0, |[ and decreases on ]|, Since this term starts 
at zero, it is positive for some time, begins to decrease at |, traverses zero at some 
point ^1 recalling that at k the value is negative, and stays negative til ^ = A; since it 

continues to decrease. Thus the sign of ^^^A^orm (C) is positive on ]0, [ and negative 

^ 2 

on ]^i,/c[, and so A^^j-j^ is increasing till and decreasing afterwards. □ 

2 2 
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Lemma 4.7. Assume k>2,alnB> and G [0, 1]. Then we have 



"^norm (0 — 



exp 



alnB 



a\nB J ^'^ 



The assumptions are already true for C = 2B when k > 4. Note that this is rather 
sharp as for ^ € [0, 1] we have X^orm (0 ^ aT' 

Proof. We use the integral representation from [Theorem 4.41 and estimate the 
polynomial hill part m'^(,^ — g) of its kernel by a simple piecewise constant function 
as indicated in the picture. We obtain 



B-e»m''{^ - g) dg 



> r B-s'' dg-m''{^-e) 
Jo 




1 



alnB 



[l - exp {-ealnB)) m^{^ - e) 



This holds for any e G [0, ^] since < ^ < ^ ensures that is increasing. So we can 
optimize e depending on ^. We obtain a suitable value when setting e by 



(4.8) 



1 — exp {—ealnB) 



k' 



To make sure that now e < ^ we use the following simple fact. 
Fact 4.9. For any i? > and r = i-<^^p(-^) the map 



[0,1}] 



1 - exp(-z), 




is bijective and increasing and for z G [0, t?] we have tz < 1 — exp(— z) < z. A 
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Let > be such that 1 - exp(-??i) = i. Namely, i?i = -In (l - i). With 
Ti = ^ '^^Pj '^i) then 1 = Ti'dik and | < ri??!. Thus we have ealni? G [O,"!?!] so that 

(4.10) nea In 5 < 1 - exp(-ea In = | < ea In B. 

k 

In particular, e < ^ follows from krialnB = ^alni? > 1. Bv lFact 4.91 with i? = ln2 



i9 

we obtain r = 2^572 ^'^'-^ ~ i) — ^^p(~t) ~ exp(— ln4) for all k > 2. Thus 
k'di = -A;ln(l -i) < ln4. Further, alnB > ^ > > This now implies 

Since ^ G [0, 1] we have the explicit expression m'^(^ — e) = ^^(^1^)1 and so 

1^= /f\ ^ ^ -fc^c ^\ _ (1 ~ g/C)^ ^ 

- toh^B^ - - alni? • • 

Since "i^i < a In B we can define t?2 by 1 — exp(— 1^2) = ^ , and according to lFact C9l 



let T2 = L_22pL_!!2) , Combining with (4.10) gives us ( 1 — I ) > ( 1 



1?2 



alnB ' 

,k / , \k 



5y — I fcriCflniJ 



exp ( — ^^^^i^Q ) and thus simplifies our above inequality to 



CAp ( ^ 

(4.11) ALrn.(0> 



alnS k\ 
By our choices 

k^x 1 



A;'i?2 > ^-T^ > 



T\T20L In a In S a In i? 

Though these inequalities get equalities with /c — )• 00, we need a precise estimate. 
Substituting /c in -A; In (l - i) < ln4 with yields 

1 / A;??! 

A;aln5-??2 = -A;aln51n 1 



T2T1 \ A;a In i? 

fcalnS, In 4 



< In ( 1 ^^^^ 1 ln4 < ln2 4 

In 4 V ka\^B . 



provided Q In 5 > □ 

Though we know the value of A^orm (^) ' is orders smaller than the above left 
lower bound. Thus let us consider A^oj-j^ on [A; — 1, A;]. 
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Lemma 4.12. If k > 3 then for ^ £ [k- l,k] we have 

■^norm (0 ^ 



alnB {k-iy. 

Proof. By definition we have Xnorm iO — Q-fcgfc+ga ■ ITheorem 4.4f s third descrip- 
tion expresses X^o,.^ on [/;: — 1 , /c] as a sum of k terms. Adding the missing term i = k 
we obtain ^t^l+^a ) noting that A is constant for ^ > fc: 

Tfc /c\ f ^ ~ p(fc-g)a cutexpfc((A:-g)alnB) 

V a In 5 J {alnBf 

To check the claimed inequality we substitute t = {k — ^)alni? G [0, alnS] (elimi- 
nating 

norm \?/ , , 

(a m is) 

We have to show that this is at least n-77— ttt which we rewrite to 



{k-l) 



0<£<fc-2 

Obviously (1 - B-^^f > (1 - e"^)'' in our situation, with equality for t = alnB or 
^ = k — 1. Thus it suffices to show for any r > 

(4.13) E 

0<e<k~2 ^ ^ ' 

The remaining proof proceeds in four steps: 

o High case: Eo<<?<A;-2 TT ^ ^• 



o Low case: > ^ii^, where i + 1 < ^i^^y^. 

o Covering: Fixing a := , ^ these cases cover all r > if A; > 4. 

y 27r(fc — 1) — 1 

o Brute-force: Prove |(4. 13)] for k = 3. (Actually, for k G {3,4,5,6,7,8,9}.) 

High case: Since e"^ > Y2o<e<k-2 Tf — ^ ^^li^ ^^se we have r > In A;. Employing 
Bernoulli's inequality (1 — T)^ > 1 — kT for T = e""^ < 1 we obtain 



(l_(l_e-)«) e--—-- <A;e 



k-l \ / ^k-1 



(A; - l)!y - V (^-1)! 

/ ^k-i \ e 



(A; - 1)! J - - ^ £\ 
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As Bernoulli's inequality is good for T close to only, it is not surprising that this 
only gives a sufficient result for r large enough. 



Low case: Assume ^ ^ > -y^^^r" where a is chosen such that ^ + 1 < -j^z^k^- 

(Since ^— ^ — is decreasing, this is always true on some interval [0,ro].) 

The condition on a implies that e"^ — + l) (k-iy. non-negative for all r > 0: 

Consider /o(t) = ^ ^t-P' ~ "I" l)- Then /g vanishes at r = /c — 1 only and thus /o 
is minimal there. The assumption on a is precisely /o(A: — 1) > 0. 

2 k 1 £ k 

Further, note that (^32)1 + (fcl)! ^ '^^ ^ ^o<e<k-i 1\ + We use this to 

obtain: 




1 



<e^-77— rTT-(i + ^)Tre^ + (i + ^ 



{k-l)l ^ ' k\ ' ' k\ {k-l)\ 



0<£<fc-2 ^ ^ 



0<e<k-2 ^ ^ ^ ^ 0<e<k-2 

— — V — — 

Notice that the value of o" only influences the set of values of r that fall in this case. 
Covering: We choose a := . = — . Recall Stirling's formula: For any n > there 

■y/27r(fc— 1) — 1 



is a ?? G ]0, 1[ such that n! = {^)'^ ^/2^^ne^^"+^ , see iRobbini (|l955l ). We thus have 



i + 1 = y^27r{k — 1) < '^cfc_^)fc-i ' a-s required. Further, we define the value Tsput where 



(fc-i) 

we split between the low and the high case: 



'21nf ifA:>8, 
Tspiit := <( 2 if 5 < < 7, 

y?-! if A: = 4. 

We start with the treatment of the cases k > 8. 

We claim that for r < Tsput we are in the low case, ie. 



,4.4, '-^^-f-if- ' 

Consider /2(t) = 1 — e""^ — i^r. It obviously vanishes at r = 0. The derivative of /2 
shows that there is exactly one maximum at r = — ln(i9), which is roughly In A;. To 
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Figure 4.3: Case coverage and Tsput 



prove (4.14) for r G ]0,rspiit] it is thus sufficient to prove that /2 is at least at the 

right boundary. First, note that by Stirhng's formula we have k\ > ^/2^^ and so 
we can estimate by |: 



'1 + a , e J 2 e 



Well, now we have 



<i 



e2 



k ■ /2 (Tsplit) = ^~ ~j: ~ k^Tspht 

>k-^-2eln-=:h{k). 

Checking that /3(e) = and the derivative f'-^{k) = (1 — is positive for k > e 
shows that fsik) > for all A; > 3. Thus for r < Tgpiit we are in the low case with 
the above choice of a. 

It remains to check that for r > Tsput we are in the high case. We use the Lagrange 
remainder estimate of the power series of the exponential function and again Stirling's 
formula to obtain 

As the left hand side is increasing in r > we consider the smallest r in question: 
let fQ{k) := f5{k,Tspiit)/k. Therein, 

(^^^y ' = exp(^-(A;-l)(^-ln2-l-lnln^+ln(fc-l) 

(I) ^ ^ 

(11) 
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Observe that the term (I) is obviously positive and increasing for k > 8. The same is 
true for the term (II): its derivative — is positive for A; > e^ ~ 7.39T and 

the value of term (II) at e^ is positive. Using this we infer that /g is increasing and 

positive. Checking /6(10) > 1 (/6(10) G ]1.19, 1.20[) now proves Eo<£<fc-2 ^^^^ > 
k for k > 10. For k = 8 and A; = 9 we just verify this inequality directly. 

It remains to consider 4 < k < 7. Here we use individual seperation positions as 
defined above: 



k 


4 5 6 7 (8) (9) 


7"split 


V7- 1 2 2 2 2 In ^ 2 In ^ 

* e e 


T 


"^split 
split 
t/l+cr 

V' fc! 


0.491 0.43JL 0.43JL 0.43JL 0.40n 0.37n 

V for low case 

0.48Y 0.40T 0.35V 0.31U 0.28U 0.25V 


0<^<fc- 


'split 

2 


4 6.33A 7.001 7.27Y 8.60JL 10.93V > k for high case 



Just check that for this Tsput the low and the high case conditions are both fulfilled. 
Summing up: the claim is proved for k > 4. 

k < 10: For A; = 3 we need an explicit check as the estimates done in the low and the 
high case are too sloppy. As the following actually is a general computational way to 



verify the inequality (4.13) we do describe it in general, show computational results 
for 3 < A; < 10 and make the critical case A; = 3 hand-checkable at the end. For the 
verification we use a small trick and brute force: First, we substitute occurrences of 
e""^ with a new variable T. The task turns into showing that the bivariate polynomial 



F,{t,T):= Y1 



(1 - (1 - T)* 



0<^<fc-2 



^fc-1 



1 



(A:-l)! 



T 



is non-negative at r = — InT for T G ]0, 1]. Now, observe that is increasing in 
r > for fixed T £ [0, 1]. If we thus replace — InT with a lower bound and we can 
show that the resulting term is still non-negative then we are done. For T E ]0, 1] 
we have — InT > 22i<i<s — T^^' This lower bound even converges to — InT, which 



actually ensures that we can always find some s that allows the following reduction. 
We consider the univariate polynomial 



9k, 



\l<^<s / 



By our reasoning, the claim follows if VT € ]0, 1] : gk.s{T) > for some s. This in 
turn is implied by 



(4.15) 



9k,si^) > A gk,s(T) has no zero for < T < 1. 



Coarse-grained integers 23 



The second statement can be checked using Sturm's theorem (ISturmlll835l ) by only 
evaluating certain rational polynomials at T = and T = 1. However, this only 
works if the chosen s is large enough. We have determined the smallest s that make 
(4.15)1 true: 



k 


3 


4 


5 


6 


7 


(8) 


(9) 


s 


4 


3 


4 


5 


6 


7 


8 




11 


13 


21 


31 


43 


57 


73 


time (sec) 


0.19 


0.29 


0.60 


2.8 


24 


235 


1704 



Though we can always divide out (1 — T)^ from g^^g the degrees are in all cases quite 
high and the computations better done by a computer. The timings refer to our own 
(non-optimized) MuPad-program used to assert (4.15) As /c = 3 is the only case 
that we do not cover otherwise we give (73^4 here: 



53,4(T)/(i-r)3 



_L(i_r) + ^r + r(i-r) 



l!i(l_T)2 + ^T2+T(l-T) 
96 ^ ^ 24 ^ ^ 

^(i_r) + — r + T(i-T) 

96 ^ ^ 24 ^ ^ 



119, 89 o 407 

1 - Tf + + T 1 - T 

144 ^ ^144 ^ ^ 

With this description you can easily see that it is positive on [0, 1[. 



n 



Finally, we put together the upper bound on A^^^.^^ . [Lemma 4^ ILemma 4.7| and 
ILemma 4. 121 in the following theorem. 



Theorem 4.16. For any k > 1 we have 

-^norm (0 ^ Cfc : = 



1 ifk = 0. 



Assume alnB > max (In 2, is^) . Then for any e € ]0, 1] with e < k - i\ or k < 2> 

2 

there is a Ck > such that for ^ £ [e,k — e] 



"^norm (0 — 



alni? 



Here we can choose 



24 Loebenberger & Niisken 



Note that is the maximum of Xnovm which in our experiments is always less then k — 

1 for k > 3. However, proving that would result in a stronger version of lLemma 4.121 
which even in the given form required quite some effort. However, we just want that 
to work for some small e and that is always granted. 

Proof. The upper bound being proven we consider the lower bound. First, keeping 
in mind that alnB > In 2 and alnB > we consider the cases A; > 3. We know 

bv lLemma 4.5l tliat X^^^^ on [e, k—e] attains its minimal value at one of the boundaries 
since G [e,k — e] (by assumption). We thus only need to consider its values at 

2 

^ = e and at ^ = /c — e. 

On the left hand side lLemma 4.71 gives 

. ^ rk ,^ / ln2 4\e'= / ln2 4\ e'^ e'' 

alnB. A_ {e) > exp ^--^J - > exp ^-^J ^ = 2" " ^ 

using the conditions on alni?. 

On the right hand side by ILemma 4.121 we find 



aln B ■ X'^^^^ {k - e) > (1 - exp (-a In S))''- — >2 



(A;-l)! - (A; - 1)! 

using again alni? > In 2. This completes the proof for the cases k > 3. 

We will now show corresponding lower bounds for the cases k = 1,2. For k = 1 
we have for e G ]0, 1] using alnB > In 2: 

alnB ■ A„of (e) = 1 — exp (—ea In i?) > 1 — exp (— e In 2) . 

Using [Fact 4.91 and e < 1, we have 

, , , l-exp(-hi2) 1 1 

1 - exp(-eln2 > P ^-e = — — e > -e. 

'-^ ' - In 2 21n2 - 2 

For k = 2 we again apply [Lemma 4.7l for the left hand side as in the general case. 
For the right hand side we show that 

a^ In-" B ■ Xl,,^ {2-e)>(^l- ^) ea In B 

for < ealnB < olni?. In 2 < olni?. Since (l — 21^) — this proves the claim. 
Now, using [Theorem 4.4! we write the left hand side minus the right hand side as 
h{ealnB,alnB) with hir,^) = exp(r - 2^9) - 2exp(r- i?) + + 1. We have 

to show that is non-negative if < r < and > ln2. The i?-derivative of /s, 

||(r,T?)=2exp(-T?)(l-exp(-T?))exp(r), 
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is positive for ■!? > 0. Thus it suffices to show that /^{t, 'd) > for the smallest allowed 
•d, which is the larger of r and In 2. If r > In 2 then we consider /^{t, t) = exp(— r) + 
2^jl^r — 1. This expression is increasing in this case (even for r > In 2 + In In 2) and 
so it is greater than or equal to /5(ln2,ln2) = 0. If otherwise < r < ln2 then we 
consider /5(t, ln2) = — |exp(T) + 215^2^ ~'~ ^ which is decreasing even for r > and 
so it is greater than or equal to /5(ln2,ln2) =0. □ 

Summing up we obtain: 

Corollary 4.17. For any e G ]0, 1] and any k > 1 we have uniformly for ^ G 
[e, k — e] 

5. Estimating the estimate A*^ 

The recurrence ILemma 3.61 for A'^' is more complex than the one for A'^ , so instead 
of solving it we estimate it. We consider also here the normed version Xnovm iO '■— 
Q,fc"gfc+ea ■ '^o better understand how the error behaves we compute it for k = 1: 

'0 if^G]-oo,0[, 



1+ga IriB J 1 \nB ,f c r- \r\ 1 r 

i+a IS-B L U if (= f1 nof 

8™ + 87ra^i+Ca lI?fc[J^,00[- 



From this we estimate Xj^oj.^ directly: 

'0 ifCG]-oo,0[, 



We have also looked at precise expressions for larger k, yet they are huge and do not 
give rise to better bounds. 

Theorem 5.1. Define values recursively by 

4 + 3 In 5,^ 
Cfc := Ck-i H 7=- (cfc_i + Cfc_i) 



for k >3 based on cn := 0, ci := (^+") ^ ^nd 

87rav B 

6 + 3a 1 4 + 3 In S _ , 

C2 := -T, 5-^= H 1=^ (ci + ci) 

87ra2 87^a^/B 

_9 + 3a 1 1 (2 + a)(4 + 31n5)hi5 

~ 87ra2 2^^a'^^/BlnB U-n'^a'^B 
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Then for any k and ^ G M we have 

Anorm iO < 4 

If a> ^= we have for k >2 and large B the inequahty 

(2'=-l)(l + a) 



Cfe < 



a 



For k = 1 the order of ci is necessarily shghtly larger. More precisely, we have for 
large B that 

[l + a) \nB 



ci < 



a 



^Jb 



Instead of defining C2 and ci wc could have left that to the recursion. But the given 
values are smaller than the ones derived from the recursion based on cq only. For 
A; = 1 the recursion would give '^"^'^^^ . For k > 2 however the improvement due to 
these explicit settings is a factor of order In 5. 

Proof. We first show that the value is bounded as claimed. For k = 1 the claim 
follows directly from the definition, since we have for B > exp(l) that {2 + a)lnB < 
2(1 + a) In B as a is positive. For k = 2 we have for In^ B < \fB the inequality 

6 + 3a 1 4 + 3 In B _ , 
l + a 1 (1 + a) In^ S 
< 3(1 + a) 



a'-^jB 



For A; > 2 we proceed inductively. We have 

^ ^ 4 + 3InS,^ ^ , 

Cfe = Cfe_l + — 7=- (Cfe_l + Cfe-l) 

OTrav-D 

^ (2^-1 - 1)(1 + a) Ing / 1 (2fe-i-l)(l + a) 

a2/B av^ Vo!lnS aln5 

^ (2fc-l)(l + a) 
a2/B 

Now we prove the remaining estimate by induction on k. The case A; = is true by 
definition of A'^ (with equality). For A: = 1 the inspection above proves the claim. The 
explicit calculation of also shows that a bound of order O / ^/b) is impossible. 
We defer the case A; = 2 to the end of the proof as most of it will be as in the general 



Coarse-grained integers 27 



case. So assume k > 3. Using the definition of A'^ from ILemma 3.6l we split A*"' into 
three summands: 

'^norm iO = / "^norm ~ S) (^Q 

Jo 



Jo 

For E(x) = ^^/xlnx we calculate as a preparative 

2E{B) InB 



(5.2) 



(5.3) 



f ^'(si+^?")lnB d^ 
^0 



aB A-Kay/B' 

4 + In B 4 + In C 



87ra\/S 87ra\/C' 

The first summand of A^^.^^ is at most c^-i by induction hypothesis. The second 
summand we estimate using (5.2) by 

InS _ _ , 

Cfc-l + Ck-l) . 



AnaVB 

The third summand is bounded by 



(ck.i + Ck-i) [' E'{B'+n^^B dg. 
Jo 



By (5.3) the third summand is at most 

4 + lnS,^ ^ , 

-= {Ck-l+Ck-l) . 

Qirav B 

This completes the proof of the case A; > 3: 

/^v InB ^ . i + lnB _ , 

\orm (0 < Cfc-l + y= (Cfc_i + Cfc_i) + y= (Cfc_i + Cfc_i, 



For k = 2 we have to improve the estimate of the first summand only, since the 
other terms anyways are of order at most O {^/^/B)■ For ^ < there is nothing to 
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prove. For ^ € [0, 1[ we find for this first summand 

Aiorm(C-^> d^< ^^^^^ ^ ^ 2 d^ 

(2 + a) In S i /■« P. 

= ^ B 2 B 2 / B 2 dg 

S-rra . L 



1-B T 



^ 2 + a 



For ^ e [1,2] we find 

Til \A ^ r (2 + a)lnS„_i_(«^ , 
X ^^""^^^"^^ ^ 8^ ^ ' ' 

JO 

(2 + a) In i3 ^ i (e-i) 
— -B 2 2 B 2 B 2 dg 



Sira 



5-1 



2 



(2-0° 

1-B 2— 



+ Xl,^{l)B-(^-'^^ f \^^dg 
Jo 



= -ii^(l-B-(C-i)°) 

a m S V / 

< — ■ — B~2 2 — -\ ' — B~2-2 -\ B~2~°' 

6 + 3a 



As A^oj^i decreases for ^ > 1 this bound also holds for > 2. Putting everything 
together the above defined value C2 bounds A^^j.^^ {^) as claimed. □ 

It is tempting to guess that we can save more Ini? factors for larger k. How- 
ever, inspecting X'^ovm shows that, say, Anoj-m (^) ^ ^ ("^)' ^^^^ ^ ^ t^' '^^ ^^'^ 

^norm (0 = 4^^^^ (v^^lns) 

6. Reestimating A'^ without Riemann 

If you do not want to assume the Ri emann hypoth esis then only weaker bounds E(x) 
on |7r(x) — Li(x)| can be used. In iFordI (l2002al lbh we found the following explicit 



bounds, the first one he attributes to a paper by Y. Cheng which we could not find. 
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Fact 6.1. o For x > 10 we have 

3/1 3 1 

|7r(x) — Li(x)| < 11.88 2;(ln x) 5 exp I (lnx)5 (Inlnx)^^ 

V 57 

o There is a constant C and a frontier xq such that for x > xq we have 
\tt{x) - U{x)\ < Cxexp (-0.2098(lnx)t(lnlnx)" 



Admittedly, these bounds only start to be meaningful at large values of x (eg. the 
first statement around 10^^^^^^). All those bounds are of the form: For all x > xq 

\tt{x) - U{x)\ < Cx{lnxy''exp{-A{lnxy^{lnlnxy'') 
' V ' 

= :E{x) 

holds. Here, C > 0, xq > 0, cq G M, ci > 0, C2 > and A > are given parameters 
(which are not always known). Note that we have that 



E{x)/ 



X 



is decreasing for large x. Actually, with the parameter sets from IFact 6.l1 this is 
already true for x > 5. Moreover, the quotient of the relative errors at x^^" and at x 

Eix'^-n'-^ (l + a)%i+°) 



EixY-^ x^E{x) 

is bounded (or even tends to zero) with x — t- cxd for any a > 0. This follows from 

^(■^)/x decreasing when a is constant, but you may also consider values for a that 
increase when x grows. 



Revisiting the proof of [Theorem 5.11 shows that only (5.2), (5.3) and the initial 
values ci and C2 depend on the specific bound E. We now use the following recursion 
for the bounds: 



Ck ■= Ck-l + 



2E{B] 
aB 



+ 1^ E'iB'+n'^^B dg^ {ck-i+ck-i 



for k > 1 based on cq = 0, and possibly values for ci and C2- 

To bound u tightly the trickiest step is bounding the integral. As our interests 

lie elsewhere we take the easy way out. We integrate by parts and use that E(^)/x 
is decreasing for the following rough estimate 

/' E'iB^^n In B d, = - + 1„ B /' d, . 

Jo aB^+" aB Jq B^+s^ ^ 

" V ' 

— B 



30 Loebenberger & Niisken 



Thus u is bounded by 



u < 




bounded 

In the following we neglect the bounded term, as we can compensate its effect for 
example by a small additional factor. Since u is small for large we expect to 
be dominated by ci = u. Precisely, for A; > 1 we have Cfc = (1 + u)ck-i + ctVaB ^ thus 

. . {l+uf-^-l k-l 



Ck = [l + uf 'u+^ ^— 1 + 



a In i? V a In 5 



Theorem 6.2. Assume that E{x) bounds \it{x) — Li(2;)| and E{x)/x is decreasing 
for X > xq and the relative error decreases fast, ie. 



{l + a)E{x^+' 



E{xy-^ x^E{x) 
is hounded under the chosen behavior of a. Then for any k > 2 and B large we have 



alni? 



a In i? 



B 



for k>2. 



□ 



This is close to optimal, we only loose a factor of order Ini? in the relative error 
compared to the used error bound in the prime number theorem: 



A: 



< 



alnB 



The assumptions on E also hold for explicit error bounds with E{x) G O [j^— ) • Due 



to the lost Ini? the res ult is only me aningful if £ > 3, so lRosser fc Schoenfeldl (jl962l ) 
does not suffice. From iDusartI ()l998l ) we can use E{x) = 2.3854^^^ for x > 355 991, 
and obtain 

_ / 1 
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7. Improvements 

We have a look at the quahty of ITheorem 3.51 when apphed to our inspiring apph- 
cation. There B = 1100 • 10^, C = 2^'^ - 1, a = ln(C)/ln(5) - 1 = 0.232T, and the 
largest k of interest is /c = 4. For these parameters we find 

1 



1 



C [0.434, 1] 



That is a great loss when we try to enclose the function of interest in a small interval. 
Actually, we can improve the theorem for the price of a slightly more complicated 
recursion. The present result was based on approximating k{g) = for p = B^^^°' G 

]S, C] by A(^) = in the recursion of IDefinition 3.ll 

fC 



B 



^B,c {0 = alnB 
We get a better bound by using 

m = 

with p = instead. 

Definition 7.1. For x > we let : 



1 



- g) ■ kiQ)B'+^'^ dg. 



+ 



InB InC 



Q, and recursively for k > 



-k-l 



1 



+ 



hiB InC 



B 



dg . 



Theorem 7.2. Write C = and fixk^f 



>0. Then for x E 

k 

k 



^>0 



we have 



,1 



17'= (x) . 



□ 



In the light of the inspiring application we now find 

k 



1 + a 



1 



C [0.957, 1] 



When we started to think about solving the recursion in IDefinition 7.11 our first trial 
was to reuse the polynomial hills m'^. Yet, that didn't want to fit nicely. Instead 
we learned from our calculations that an exponential density instead of a linear one 
would be easier to connect to the polynomial hills. So we tried to approximate like 
this and got 



Inp {1 - g)lnB + gliiC 
The exponent in fj is chosen such that for g = and ^> = 1 we have equality. It turns 
out that this approximation is even better than the one before and at the same time 
easier to handle. Thus we replace the functions T^^ with another family rj^: 
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Definition 7.3. For ^eRweletrf := ^, and recursively for k > 



iO := a In S / - g) B^+^" dg . 




Theorem 7.4. Write C = B^+°' and Ex k £ N>o. Then for x G M>o we have 

k 



ln(l + a) 



a 



(1 + a)Mi+^)' 



1 



rj'' {x) 



Proof. To prove this we have to relate the function k: [0,1] — t- M>o, q i— >• 
1/ln {B^+S'^) occurring in the definition of to the function i) : [0, 1] — y IK>07 Q ^ 
fj{Q) replacing it in the definition of rj''. Routine calculus shows that the function 
k/fi is at most 1, namely at g = and g = 1, and assumes its minimum value 



ln(l+a) 
a 



(1 + a)i"(i+") at g 



a— ln(l+a) 
a+ln(l+o) • 



□ 



Testing this with the parameters a = 0.232T and k = 4 from our inspiring application 
we obtain 



ln(l + a) 



a 



(1 + q)Mi+") \ I 



C [0.978,1]. 



To get a better impression we have plotted the lower interval boundary as a function 
of a for all three cases in Figure TTT] We see that for small values of a we obtain 
good approximations of ^ and all our attempts give only weak results for large a, 
but the one with fj is always best. 

If we now rewrite the recursion to one for 



(0 



, , /-f ln(l + a)^ 



we find that rj^ 



-^^norm- ^o We'll obtain the solution from the polynomial hills 



as in [Theorem 4.41 for A'^: 



^norm (0 ' ^ 



-g{alnB-\n(l+a)) ~k 



m^{i - g) dg. 



The only difference is that instead of a we have a — — j^-^- With this replacement 
[Theorem 4.41 becomes: 
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Figure 7.1: Comparing the quality of A, V, rj: the top pictures show the integral 
kernels and their ratios for a specific a, the lower pictures show the minimum of the 
ratios as a function of a 



Theorem 7.5. For any ^ G M we have 



Jo 







5 / 

-—] m\i-Q)dQ 
1 + a J 



CI InC 
' B/lnB 
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'/norm 



(-alnB + ln(l + a))*^ 



^ (-aln5 + ln(l + a))^P|-^-im^(0|, 

0<^<fe-l 



~k AV .Nj Cutexpfc j-j^-i) {alnB - ln(l + a))) 

'/norm \S/ ^ I / ^ / 1 D i Wi i ' 

o<i<L€J (-alnB + ln(l + a)) 

wiere cutexpjt(C) = exp(C) - Eo<£<fe-i S = E£>fc 

8. Non-squarefree numbers are negligible 

Considering ^ (x) we immediately observe that the ordering of the counted prime 

lists are not important and we can group together many such elements. To get a 
precise picture, we define the sorting of a tuple P = (pi, . . . ,pk) in the following way: 



S{P) ■.= (^{j<k \ rankj P = i} 



i<k 



where rankj P = # {p( \ £ < k A p( < pj}. Now the count ^ (x) can be partitioned 
using the sets 

As{x) := |P = ip^, . . . ,p,) G (P n ]B, C]f I ^ = j'^^ 'f' - , 

namely k% (x) = Es^^s(^) where S runs over all possible sortings. The above 
intuition would imply that many of these sets are essentially equal. We group them 
by their type 

T{S) := (#5,),<, . 

Given any type T = (Ti, . . . ,Tr), there are exactly (^) = rp^i 'j, , different sortings S 
of type T. This corresponds to possible reorderings of a specific vector P G As{x) 
for a sorting S of type T. The type of such a vector P is defined to be T(S). It is 
clear that the type of P is invariant under permutations, yet not its sorting. 

Lemma 8.1. Let T be a type for k elements. 

(i) There exists a sorting S{T) of type T such that all vectors in As(^t){x) are 
increasing. 
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(a) IfT{S) = T then there is a permutation a of k elements such that for all x we 
have As{x) = As(^T){xT . 

(Hi) More precisely, for any sorting S of k elements the following are equivalent: 

(a) T{S) = T. 

(b) 3a: S = S{Ty. 

(c) 3(7 : Vx: As{x) = As^r^i^T ■ □ 
Noting that # {5 | T{S) = T} = we have 

4,c (^) = E E #^^(^) = E (t) *Mt){x) 

T S:T{S)=T rjp \ / 

On the other hand we have i:^ ^ (x) = ^ ' H^^S(T){x)- In particular, we can 
deduce 

Actually, for large B (and C and x) we have 7r|^ (j (x) ~ ■^^s c ('^)- This stems from 
the following fact that ^As{x) is asymptotically much smaller than #^s(i,,,,,i)(x) 
for any sorting S of k elements of type different from (1, . . . , 1). 

Lemma 8.2. For any sorting S of k elements of type different from (1, . . . , 1) there 
is a sorting S' of k — 1 elements such that we have ^As{x) < i^As'{x/ B) < 



x_ 
B- 



Proof. Take S as specified. Let i be a position which does not occur as a singleton 
in S. Further, say t £ Sr, and let r = > 2. Let S~ be the sorting with Sr 

removed, and S' the sorting with t removed. (Retaining the old indexing is easier, 
yet then indices run over {1, . . . ,k} \ Sr, or {1, . . . , A:} \ {t}, respectively.) Then 

#As{x)= E #A',^{x/pI)< E #A',4x/ptB) = #As'{x/B). 

Pt&r\]B,c] pt&n]B,c] 

Here, A' denotes a variant of A consisting of prime vectors where at the positions j 
which should have smaller primes than pt originally still satisfy that, and same for 
positions with primes larger than pt. For the inequality note that S~ = {S')~ . Since 
obviously ^As{z) < z we are done. □ 

Combining this with ^As{x) = c (•^) ~ '^B C (•^) ^ ® (iikb) ' shows that there 
must be a large summand, which can be only #^5(1 ;^)(x). 

The number s{k) = (^) of sortings of k elements is called ordered Bell number. 
We can also recursively define them: s(0) = 1, s{k) = X^o<r<A:-i (r)^('^)- According 
to IWil j (Il994h . page 175f, we have s{k) = ^^^+1^ + O ((0.16)'= A;!) . In particular, 
s{k) is small in comparison to 2^~^k\. Using [Lemma 8.2] for a comparison yields the 
— now immediate — following 
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Lemma 8.3. We have 



Proof. 



k\ ■ 7r^_c (x) - (x)| < Et - (t)) f = (^!2^'' " ^(^)) §• ° 



Compared to the error bound in q {x) — c' i-^) — c (•^) ^ ^("^^ ^^^^ 
negligible when B is large. Here we assume k > 2 since the present observations are 
irrelevant for /c = 1, namely tt^ q = k}^ q- 



9. Results on coarse-grained integers 

We are e^oing to combine the results of the last section with [Theorem 4.16) and The- 
orem 5.1 to finally arrive at our main theorem. 

Analogously to IDefinition 3.H we define an approximation function for the func- 
tion vr^ ^ (x). 

Definition 9.1. For x > and k>Q we define 

^B,c (x) ■■= l^^'h,C (^) aiid (x) := ^K%^c (^) + 2'="'^- 
Similarly to IDefinition 3.11 we can also recursively define tt^ ^ (x) by 



if X < 1, 

1 if 1< X, 



T^B,C (x) 



k 



c4;^(xM) 



It is also possible to define vf^ ^ (x) similarly based on IDefinition 3.11 We can now 
describe the behavior of vr^ ^ nicely and give our main result. 

Theorem 9.2. Given x G M>o and /c G N. Then the inequality 

T^B,c{x)-^%^c{x)\ <^%^c{x) 

holds. 

Proof. Bv lLemma 8.31 we have 



vr' 



B,C (^) ~ '^'^■B.C 



X 



Thus using the triangle inequality and lTheorem 3.2] we obtain 



A,C (x) - l^'^B,C (x) 
which proves the claim. 



< ^^B,c (x) + 2 



fc-i 



B' 



□ 
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Theorem 9.3. Fix k > 2. Then for any e > and B tending to infinity, there are 



for X e [^''(1 + e),C''(l -e)] values s,se 



i 1 

(l+a)* ' 



such that 



Proof. Bv IDefinition 9.11 we have 



ITheorem 3.51 tells us that there is a value s E [i+a)'' ' ^'^'^^ that 

implying that for the same value s we have 

^B,C i^) = ^A*^ (^) • 
Considering vf^ (x) we have by IDefinition 9.11 that 

^B,c (^) = ^4,c {x) + 2''-'^. 



Applying [Theorem 3.51 gives a value s G 



1 



which directly gives 



This proves the theorem. 



(l+a)*^ ' 



B 

such that 



k\ 



B 



n 



Unrolling our results on A'^ {x) and A'^ (x), namely [Theorem 4.161 and ITheorem 
gives a slightly weaker result. 

Theorem 9A. Let B < C = B^+" with a > and fix k > 2. Then for any 
(small) e > and B tending to infinity we have for x G [i?'^(l +e),C^{l — e)] a 



value a G 



fc!(l-|-a)'= ' k\ 
.k 



with Ck = min ( 2 



-As'' n-fc £ 



fc-2 



sucil tiiat 



X f^t- — 1 ^ 

-^ + 2 "s- 

B B 
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Proof. Bv lTheorem 9.3l we have values s.'sG 



(l+a)fe 



such that 



TT' 



B,C 



kl 



B 



Bv lTheorem 4.161 we have for e small enough that 

A'=(x)G [4,1] 

and bv [Theorem "Othat 



a X 
InB 



This gives the claim. 



10. Numeric evaluation 



□ 



To discuss the quality of our results we consider again the example parameters B = 
1100 • 10^, C = 2^7 - 1, a = ln(C)/ln(5) - 1 = 0.232T, k = 4 from our inspiring 
application. For = 2, A; = 3 we can give similar pictures; of course, the errors 
are even smaller in these cases. At present we do not have efficient algorithms for 
computing q itself. However, based on our estimates we can compute values for 
encapsulating intervals in three variants, listed in increasing quality: 

o A estimate (jTheorem 3.5P : 



1 



(1 + 



r] estimate (jTheorem 7A\i : 
'ln(l + a) 



-i-^(l + a)i"(i+o,) (x) -Cka^x, rf{x) + Cka^x 

a J 



o n 



estimate (jTheorem 3.2p : 



4,C (x) - K%n (x) , 4,C (x) + (x) 



The K estimate was easiest to obtain and is of course the most accurate one, however, 
it is difficult to evaluate. The A estimate was easy to obtain and compute. But it is 
of course the least accurate of the three. The r] estimate was slightly more difficult 
to find, is as easy to evaluate as the prior, and it is much more accurate. As usual 
we write x = B^~^^°^ and use ^ as a running parameter. 
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0.2 • 10 



0.1 • 10 



0-10 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 ^ 



Figure 10.1: Absolute behavior of the estimates for {£,) ■ The hght gray area 
shows the A estimate, the dark gray area the r] estimate, and the black area (weh, 
yes) shows the k estimate. The parameters are B = 1100 • 10^, C = 2^^ — 1, a = 
ln(C)/ln(5) - 1 = 0.232T, k = 4. 



0.3 • 10- 



0.2 • 10- 



0.1 • 10- 



0-10- 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 ^ 

Figure 10.2: Behavior of the estimates relative to x = B^^^°'. Colors and parameters 



are as in Figure 10.1 



Figure 10.1 shows the absolute behavior of all estimates. We observe that the 



absolute errors at the right margin are huge. This is expected as also the error 
estimates in the prime number theorem only bound the relative error. However, the 
picture completely conceals information about the middle and the left part of the 
interval [0, k]. 

To see more we divide hy x = B^^^'^ and therefore obtain estimates for the ratio 



of C] -grained integers x in Figure 10.2 This reveals a lot about the quality of 
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our estimates. The black area indicates the best that we could hope for, namely the 
estimate based merely on lPrime number theorem 0.31 However, as this is difficult to 
evaluate we have to approximate once more. The A estimate, shown in light gray, is 
clearly only of use to get a rough idea. The r] estimate, however, is rather close to the 
actual behavior and may well serve as a basis for stochastic fine tuning of algorithms 
like the general number field sieve. 

Last, Figure 10.3 illustrates the size of the various errors terms relative to k^: 



A error: (l - j^Taj^) (^)- 

terror: ^1 - (^ '"^^+"^ (1 + a)T^^I^~°)^^ (x). 

A error bound: Cka^x. Unconditional A error bound: c^a^x. 

A error: A'^ (x). Unconditional A error: A^ (x). 

K error: K^f^{x). Unconditional k error: k^(^(x). 

Non-squarefree error bound: (2''^^kl — ^. 
Non-squarefree error: (2'^~^A;! — s(A;)) 



For the unconditional errors we use iDusartl 's unconditional bound on |vr(x) — Li(x)| 
which is given by E{x) = 2.3854^ for x > 355 991. 

The figure shows that all the error terms but the A error are sufficiently small. 
Our best choice is the rj estimate which is ruled by the 77 error and the A error. 
Both are fairly less than 3% of the target value at least in the middle of the interval. 
The estimations are more difficult close to the boundaries. It is also positive that, 
at the parameters of our interest, most error terms are still comparative in size to 
the contributions of the /? error, which is induced by the prime number theorem. 
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Definitely, tlie rj estimate, combining the 7] error and the A error, is good enough for 
practical purposes, as for example the fine tuning of the general number field sieve. 
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